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ABSTRACT 

We have studied the OH masers in the star forming region, W3(OH), with data obtained from the Very Long 
Baseline Array (VLBA). The data provide an angular resolution of ^5 mas, and a velocity resolution of 106 
m s" 1 . A novel analysis procedure allows us to differentiate between broadband temporal intensity fluctuations 
introduced by instrumental gain variations plus interstellar diffractive scintillation, and intrinsic narrowband 
variations. Based on this 12.5 hours observation, we are sensitive to variations with time scales of minutes 
to hours. We find statistically significant intrinsic variations with time scales of ^15-20 minutes or slower, 
based on the velocity-resolved fluctuation spectra. These variations are seen predominantly towards the line 
shoulders. The peak of the line profile shows little variation, suggesting that they perhaps exhibit saturated 
emission. The associated modulation index of the observed fluctuation varies from statistically insignificant 
values at the line center to about unity away from the line center. Based on light-travel-time considerations, the 
20-minute time scale of intrinsic fluctuations translates to a spatial dimension of ^2-3 AU along the sight-lines. 
On the other hand, the transverse dimension of the sources, estimated from their observed angular sizes of about 
^3 mas, is about 6 AU. We argue that these source sizes are intrinsic, and are not affected by interstellar scatter 
broadening. The implied peak brightness temperature of the 1612/1720 maser sources is about - 2 x 10 13 K, 
and a factor of about five higher for the 1665 line. 

Subject headings: interstellar: molecules - masers - radiation mechanism: non-thermal 



1. INTRODUCTION 

Interstellar hydroxyl maser sources are found in the Galaxy 
and many external galaxies. In the Galaxy, these sources are 
associated with star forming regions (predominantly in the 
1665 and 1667 MHz lines), IR and late type stars (mainly 
at 1612 MHz) and supernova remnants (exclusively at 1720 
MHz). 

Numerous VLBI observations (e.g., Gwinn et al 1988; Fish 
et al 2005; 2006) have indicated that the apparent angular 
sizes of OH masers increase with distance in the Galaxy. OH 
maser sources in the inner Galaxy show larger angular sizes 
compared to sources in the outer Galaxy. These facts strongly 
suggest that the observed angular size is strongly dependent 
on interstellar scatter broadening. Some years ago, Burke et 
al (1968) and Gwinn et al (1988) have proposed that the ap- 
parent broadening of masers shows a wavelength dependence 
oc A 2 , based on OH maser sizes at 1665 MHz and H2O maser 
sizes at 23 GHz. These data suggest that interstellar scattering 
is the dominant cause for the observed angular sizes of these 
sources. 

Desai, Gwinn & Diamond (1994) used the VLBA to study 
the details of the interstellar broadening of the OH masers in 
the distant (~14 kpc) HII region W49N. Anisotropic broaden- 
ing was observed with the orientation of the minor axis prefer- 
entially parallel to the Galactic plane, an additional confirma- 
tion of the influence of interstellar scattering on the scattered 
sizes. An unambiguous determination of the intrinsic sizes 



of OH maser sources would represent an important input in 
understanding the OH maser pump and emission mechanism. 

In addition to the intrinsic sizes, an important property 
of the masers that provides information about the pumping 
mechanism is the variability of the maser intensity as a func- 
tion of time. Variability on long time scales (weeks to months) 
has been discussed by numerous investigators (Schwartz, Har- 
vey & Barrett 1974; Coles, Rumsey & Welch 1968; Zucker- 
man et al 1972; Rickard, Zuckerman & Palmer 1975; Gru- 
ber & de Jager 1976; Clegg & Cordes 1991). The sugges- 
tion was made that the observed variability may arise from 
changes in the number density of molecules exhibiting maser 
phenomenon, or physical conditions in the maser column. In 
contrast to the long timescale variations, short timescale varia- 
tions have also been detected (Zuckerman et al 1972; Rickard 
et al 1975). Salem & Middleton (1978) provide a model con- 
sisting of a sudden onset of a pumping mechanism that could 
cause rapid quasi-periodic fluctuations in the observed inten- 
sity. The predicted fluctuations on time scales of roughly a 
day would have a ^25% modulation index. 

Evans et al (1972) have investigated eight well known OH 
masers with the goal of describing the statistics of the radi- 
ation. They sampled the output of the 140 feet NRAO tele- 
scope rapidly and find that the radiation is of Gaussian nature 
to within a level of about one per cent. 

The major existing study to date of short term variations 
consists of Arecibo (beam size of 3 arcmin) and VLA (beam 
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size of 15 x 5 arcsec) observations probing variability on time 
scales in the range 16 s to 2 hours (Clegg & Cordes 1991). 
Typical variations were detected at the 5-10 per cent level 
with some large variations at the 100 per cent (or more) level 
for the sources W75S and NGC 6334F. As Clegg & Cordes 
point out, an identification of any intrinsic variability in these 
sources would provide a unique opportunity to determine the 
source extent using light-travel-time arguments. They do find 
prominent variations in some of the sources with time scales 
of ~20 minutes (e.g S269). However, it is very difficult to 
distinguish between variations that are intrinsic to the source, 
and those that might arise from interstellar diffractive scintil- 
lations. These authors do not decisively choose between an in- 
trinsic or scattering origin for the fluctuations. However, they 
do stress that if the fluctuations are due to interstellar scintilla- 
tion, implausible brightness temperatures would be required. 

As pointed out by Clegg & Cordes, several aspects can af- 
fect their conclusions: 

• The 3 arcmin beam size of Arecibo, or even the 15x5 
arcsec beam of VLA imply that a number of individual 
maser sources are observed simultaneously. Thus the 
observations can consist of incoherent superposition of 
multiple sources. 

• The velocity resolution of the VLA was only 1.1 km 
s" 1 ; thus in some cases multiple velocity components 
may exist in a single velocity channel since a spectral 
resolution at the level of < 0.1 km s" 1 is required to re- 
solve complex OH maser lines. 

• The total time span of each observation was only ~ 2 
hours, thus limiting the information concerning inten- 
sity variations to time scales of about an hour. 

In the current study, we have circumvented these problems 
by using data from the Very Long Baseline Array (VLBA) 
of the National Radio Astronomy Observatory 1 to produce 
time series of OH observations at all four spectral lines of 
OH in the source W3(OH). With a total angular extent of the 
maser region of ~ 2 arcsec, many maser spots are observed 
simultaneously. With the ~ 5 mas spatial resolution and a 
spectral resolution of ~ 0.1 km s" 1 , the confusion problem is 
solved since strong isolated maser sources can be observed at 
unique positions and velocities in three of the four OH lines. 
With the snap-shot capability of the VLBA, images are made 
at intervals of 1 min over a 12.5 hour period. The modu- 
lation indices of the maser lines are derived as well as the 
power spectra. After summarizing the observations in ^2] our 
analysis procedure is described in |j3] In |4] we discuss the 
possible contributors to apparent intensity variations, and de- 
scribe our model incorporating their spectral characteristics. 
We have adopted two methods to explore the intrinsic vari- 
abilities in these maser sources; these are summarized in ^5] 
and 5j6] As we describe in these sections, we find signifi- 
cant intrinsic variations. Implications of these variations are 
described in the last section, Throughout this paper we 
assume that the distance of W3(OH) is 2.0 kpc based on the 
recent VLBA parallax determinations of Xu et al (2006) and 
Hachisuka et al (2004). 

1 The National Radio Astronomy Observatory (NRAO) is a facility of the 
National Science Foundation operated under a cooperative agreement by As- 
sociated Universities, Inc. 
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FIG. 1. — The total u-v coverage of the data set at 1665 MHz. The two 
axes are in units of kilo-wavelengths. 



2. OBSERVATIONS 

The dataset is the observations analyzed by Wright, Gray 
& Diamond (2004a,b and 2005; hereafter WGDa,b,c). They 
observed W3(OH) on 1996 August 2 using the VLBA in all 
the four ground-state lines (1612, 1665, 1667, & 1720 MHz), 
simultaneously. We have obtained the calibrated u-v data from 
Phil Diamond and analyzed the data in order to search for 
short term time variations. The data were recorded with full 
polarization information with a bandwidth of 62.5 kHz at all 
four OH maser lines. The assumed rest frequencies of these 
four spectral lines were 1612.231, 1665.402, 1667.359, and 
1720.530 MHz, respectively. With 128 spectral channels, the 
channel separation is 488 Hz, and the resolution of 586 Hz 
corresponds to a velocity resolution of 0.11 km s -1 at 1665 
MHz. 

A full description of the adopted calibration procedure can 
be found in WGDa. For the purpose of calibration, the sources 
3C84 and J 16 11+343 were observed during the run. In par- 
ticular, the amplitude calibration was carried out using the 
VLBA parameters (gain curves and system temperature de- 
terminations) known a priori, and in addition using the "tem- 
plate fitting" method using the AIPS task ACFIT. The overall 
flux density is thus determined with an accuracy of a few per- 
cent. The phase calibration was carried out using a phase ref- 
erence velocity channel at 1665 MHz at the velocity of -47.46 
km s -1 from W3(OH). This calibration (at the sub-minute time 
scale) was carried out using FRING in AIPS followed by a 
self-calibration in order to remove any effect of source struc- 
ture. No amplitude self-calibration was performed. The phase 
corrections were then applied to all channels. In Figure 1 the 
total u-v coverage of the data is shown for the data at 1665 
MHz. 

3. ANALYSIS PROCEDURE 

The first step in the data analysis is to identify bright 
sources that are spatially well isolated. At the end of this step, 
three sources were chosen, whose properties are summarized 
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in Table 1 . The spectra from each of the three positions are 
shown in Figure 4 for the 1612, 1665 and 1720 MHz lines 
based on the 12.5 hour observation. At 1667 MHz no strong 
isolated source could be identified that did not exhibit a pro- 
nounced velocity gradient in adjacent channels. We have not 
analyzed the 1667 MHz data since the gradients may be the 
result of blended sources. The three positions and six lines 
(RR and LL for each) shown in Table 1 were chosen for fur- 
ther analysis. 
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02:27:03.818 


+61:52:24.439 
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1612 


02:27:03.818 


+61:52:24.439 


RR 
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1665 


02:27:03.825 


+61:52:24.653 


LL 




1665 


02:27:03.825 


+61:52:24.653 


RR 


3 


1720 


02:27:03.829 


+61:52:24.704 


LL 




1720 


02:27:03.829 


+61:52:24.704 


RR 



Table 1 

List of maser spots used for the present analysis. The 
second column gives the line tag in mhz, the third 

and the fourth give the right ascension and 
declination (j2000) of the sources, along with the 
circular polarization tag indicated in the last 

COLUMN. 

The positions for the six lines were obtained after correc- 
tion of the WGDa,b positions. Since the August 1996 data is 
not phase referenced and is self calibrated, the absolute co- 
ordinates are uncertain at the 0.1 arcsec level. The absolute 
positions were tied by WGD a,b to earlier determinations of 
absolute positions by Gray et al (2001) with an uncertainty of 
— 10 mas. We have used the corrections to the WGDa,b po- 
sitions published by WGDc. In addition, in December 2004, 
Palmer and Goss (unpublished) have carried out phase ref- 
erenced VLBA observations with an astrometric precision of 
~1 mas. These authors have determined that the corrections 
to the WGDa,b coordinates are +48 ± 5 mas in right ascen- 
sion {i.e. add 48 mas of arc to the coordinates published by 
WGDa,b) and +184 ± 5 mas in declination. WGDc have sug- 
gested that these values are +54 mas and +175 mas, respec- 
tively. Both determinations are in reasonable agreement. The 
coordinates in Table 1 are the WGDa,b values using the cor- 
rections determined by Palmer and Goss. 

Numerous trial snap-shot images were made in order to de- 
termine the minimum time interval over which a successful 
and reliable image can be obtained. For these isolated strong 
sources (-10 Jy beam" 1 ), we find 1 minute as the shortest in- 
terval at which reliable images can be constructed with min- 
imal problems of confusion (see below for a special case for 
the 1612 MHz line). In the top panel of Figure |2] we show 
the u-v coverage of a one minute interval at an hour angle of 
4 hours for the 1665 MHz observation. To give an idea of 
the quality of map, we also show the map produced by a one 
minute observation in the bottom panel. Moreover, in the top 
and the bottom panels of Figure|5] we show the cross section 
of the beam along the right ascension and the declination di- 
rections, respectively. In order to demonstrate the signal to 
noise ratio of the detection of the spectral lines, we also show 
in Figure0]the spectrum of the 1665 MHz line corresponding 
to the 1-min integration. 

Out of the three sources given in Table 1, the sources at 
1612 and 1720 MHz have been identified as Zeeman pairs. 
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FIG. 2. — Top: Typical u-v distribution with 1-min integration, correspond- 
ing to the 1665 MHz observation at an hour angle of 4 hours. Bottom: Typical 
one minute map produced with the 1665 MHz data. 



The corresponding longitudinal magnetic field strengths are 
given in the table from WGDa. At 1665 MHz, we do not ob- 
serve a statistically significant shift between the line profiles 
in RR and LL channels. The integrated (over 12.5 hours) line 
profiles of these Zeeman components are shown in Figure|5] 

The data analysis to identify and compute the line profiles 
was carried out using MIRIAD. Since some of the sources are 
clearly Zeeman pairs with velocity separations of the order of 
— 1 km s" 1 , (several freq channels), we carried out the anal- 
ysis in both RR and LL polarization channels independently. 
The rationale for this procedure is that the intensity variations 
observed in the two polarizations may not necessarily be cor- 
related even for small (< 0.1 to 0.2 km s" 1 ) velocity separa- 
tions. The full 12.5 hour observation was used to constrain the 
source position and shape parameters (major and minor axis 
and PA). Then, we determined the flux density at 1 minute 
intervals, with the source position and the shape parameters 
held fixed. For this purpose, the MIRIAD task UVFIT was 
used. At the end of this step, we obtained line profiles ev- 
ery minute in RR and LL channels separately for all sources 
listed in Table 1. These profiles were arranged chronologi- 
cally, forming a dynamic spectrum of intensity as a function 
of time and frequency (or velocity), for each source, in RR 
and LL separately. 

For the 1612 MHz line, it was necessary to carry out a mi- 
nor correction due to source confusion. The major 1612 MHz 
line has flux densities of 15 Jy (RR) and 4 Jy (LL). There is 
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Plot file version 5 created 26-SEP-2005 12:53:50 
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FIG . 3 . — Top: Beam pattern as a function of right ascension and declination 
are given in the top and the bottom panel, respectively, for a data of one 
minute duration. 



another source nearby (displaced by -12 mas in right ascen- 
sion and +1 mas in declination) in the same velocity range. 
The flux density of this source is ~1 Jy, which is about 10% 
of the more intense source. This weaker confusing source was 
imaged using the 12.5 hour observation; the mean flux density 
was subtracted from the entire u-v data base before the time 
series was constructed for the two major lines at RR and LL 
polarization. This method is based on the assumption that the 
time variation of the confusing source is minor compared to 
that of the brighter source based on the relative weakness of 
the confusing source. 

4. ORIGIN OF FLUX DENSITY VARIATIONS, AND OUR MODEL 

Flux density variations with time scales of tens of minutes 
can be introduced by various effects that are independent of 
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FIG. 4. — Spectral line profile of W3(OH) corresponding to the one minute 
integration at 1665 MHz. 



intrinsic variability. In order to explore the possibility that 
these maser sources exhibit intrinsic variability, it is essen- 
tial to understand the nature of additional extrinsic sources of 
variability on relevant time scales. 

First, any uncorrected instrumental gain variations in the 
receiver system will introduce apparent intensity variations. 
However, these changes will be correlated across spectral 
channels. In other words, we expect these changes to have 
a correlation bandwidth that is far wider than the typical spec- 
tral line widths of the maser sources of ^0.5 km s" 1 . 

There is an effect that can potentially introduce flux density 
variability on short time scales. This effect arises from inter- 
stellar diffractive scintillations; two of the important parame- 
ters characterizing it are, namely, the decorrelation bandwidth 
(At/), and the diffractive scintillation time scale (T^f). The lat- 
ter ranges typically from minutes to hours. If the decorrelation 
bandwidth value is comparable to, or narrower than, the spec- 
tral line width of the maser source, i.e. less than a few channel 
widths or less than 0.2 to 0.4 km s— 1, then we would consider 
the temporal variations introduced by diffractive scintillations 
to be "narrowband", and would expect significant differential 
variation within the line profile. In order to assess the sit- 
uation correctly, we must estimate the expected value of the 
decorrelation bandwidth of diffractive scintillations along this 
sight-line. 

Since these maser sources are spectral line emitters, mea- 
surements of the decorrelation bandwidth are often uncertain. 
However, the measured angular diameters of these sources 
provide an indirect estimate. The apparent angular width of 
many of the sources in W3(OH) are measurable in this data 
set, and is typically ^3 mas (Palmer & Goss, Private commu- 
nication). If we assume that this width is exclusively domi- 
nated by interstellar scatter broadening, i.e. the intrinsic an- 
gular size is significantly smaller than the apparent size, it 
is possible to calculate the expected relative time delay as- 
sociated with the scattered rays with respect to the direct path 
(Gwinn et al. 1993; Deshpande & Ramachandran 1998). This 
effective time delay is identical to the characteristic temporal 
broadening scale (r sc ) in the case of pulsar pulse profiles. At a 
distance of 2.0 kpc, and assuming that the scattering material 
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FIG. 5. — Integrated line profiles of W3(OH)corresponding to the three sources summarized in Table 1. Profiles with solid lines correspond to LL polarization, 
and the ones with 'dashed' lines correspond to the RR polarization. 



is uniformly distributed along the line of sight, the delay is 
16 c lnz 

where D is the distance to the object, 6 H is the angular width 
of the source (full-width at half maximum), and c is the speed 
of light. Since the decorrelation bandwidth and the effective 
time delay obey the following "reciprocal" relation, 

27TT SC 5v w 1, (2) 

the value of the decorrelation bandwidth, 5v is then predicted 
to be ~35 kHz. This is significantly greater than the typical 
line widths (1-2 kHz, or 0.2-0.4 km s" 1 ). Of course, this es- 
timate is a worst-case estimate, where we have assumed that 
the observed width entirely arises from interstellar scattering. 
Even with this worst-case estimate, the expected differential 
variation within the spectral line profile as a function of time 
is only about ^1% or less. 

Apart from the above mentioned causes for intensity fluc- 
tuations, any additional observed variations may be assumed 
to be intrinsic to the maser source. Of course, we have no a 
priori knowledge about the nature of the intrinsic variations 
in OH masers. However, for the current analysis, we assume 
that the correlation bandwidth of any intrinsic variations is 
small enough to consider them as "narrowband". That is, 
such variations in one channel would, in general, be uncorre- 
cted with those in the other channels separated by our veloc- 
ity/spectral resolution or more. There is an important reason 
for making this assumption. The possible intensity variations 
due to instrumental gain instability and interstellar diffrac- 
tive scintillations are expected to be correlated across a veloc- 
ity/spectral range much wider than the maser line widths, and 
hence can be treated as "broadband" variations, i.e. as mod- 
ulations that are common to all the observed spectral/velocity 



channels. Therefore, in our analysis, we will be sensitive to 
only those intrinsic variations that are narrowband in nature, 
thereby clearly distinguishable from the instrumental and in- 
terstellar effects. 

In order to quantify the observed variations, we have 
adopted the following procedure. We make a distinction be- 
tween the different contributors to the observed flux density 
variance, 

^obs(v) = a 2 h + af(v) + a 2 n (v). (3) 

The three terms in the above equation result from a "broad- 
band" modulation, possibly intrinsic narrowband variability, 
and measurement uncertainty, respectively. For this purpose, 
we model the observed intensity S o/K (f,v), i.e. a function of 
both time (f) and velocity (v) as 

SoUt , V) = S aV e(v) [1 + f„(t)] [1 +/ v (v,f)] + «(v,f), (4) 

where fb and f, are the zero-mean fractional variations that 
are of "broadband" and "narrowband" nature, respectively, 
and n is the measurement noise, assumed to have zero mean. 
The function 5 ave (v) refers to the average line-profile as a 
function of velocity. The time averaged cross-correlations 
between f„ f, and n are expected to be zero, i.e. they are 
assumed to be mutually uncorrected. The parameter fb, the 
"broadband" modulation, is independent of velocity (or fre- 
quency), and is therefore described as a function of time 
alone. Adopting the above formulation, the observed variance 
a 2 bs in a given velocity channel can be expressed as 

^obs(v) = Si e (v) [a 2 c + <7 v 2 (v)] + <t„ 2 (v), (5) 

where <r c and <7 v (v) are the root-mean-square (RMS) fluctua- 
tions characterizing the fractional temporal variations fb and 
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/ v (v), respectively. a„(v) is the RMS uncertainty in the mea- 
surements, which can be estimated in a straight-forward man- 
ner as -\/^e 2 (v,f)/^V or \f (e 2 {v,t)),. Here, (x), represents 
the average of "jc" over time, and e(v,t) (l-a error), is avail- 
able for each of the N (typically, 300-400) time samples in the 
dynamic spectrum S b s (t,v) at a given velocity v. The e(v,t) 
estimate for relevant channels also includes the contribution 
of the maser emission to the system noise. The desired quan- 
tity <j v (v), the modulation index associated with the "narrow- 
band" variation (uncorrelated across velocity channels), can 
be estimated if a c is known. 

5. CORRELATION ANALYSIS TO ESTIMATE "BROADBAND" 
VARIATION 

In order to estimate and remove the possible contribution 
of any broadband variation (characterized by a c ) from the ob- 
served variance, we examine the cross-correlations between 
the intensity fluctuations in all velocity-channel pairs. Auto- 
correlations are excluded since they, as seen from equation 
5, are contaminated by the variance of measurement noise 
and narrowband fluctuations. Since fy, f v and n are mutu- 
ally uncorrelated, the cross-correlation, expressed as "cross- 
variance", between intensity fluctuations in any pair of chan- 
nels (about their respective mean intensities) is given by 

CobsOl ) v l) = {Sdtaif, Vl ) S obs (f , V 2 )), 

= Save(Vi) S ave (v 2 ) ((MO fb(t))t + 

(/v(vi,0/v(v2,0>t) + (n(v u t)n(v 2 ,t)) t 
= S ave (vi) 5 ave (v 2 ) (f b (t) f b (t))t 

= SaveCVl) Save(v'2) of> (6) 

where vi v%. In a given velocity channel, the contribu- 
tion due to any "broadband" intensity fluctuations, treated as 
an amplitude modulation, is expected to be proportional to 
the mean intensity associated with that channel. The cross- 
correlations between such fluctuations then will be propor- 
tional to the product of mean intensities of the corresponding 
pair of channels. The possible contribution from any "nar- 
rowband" fluctuation, f v (v,t), will also share this proportion- 
ality, but the associated cross-correlation is expected to have 
a zero-mean value. Similarly, the cross-variance of the mea- 
surement noise in any two different (and independent) veloc- 
ity channels is also expected to be zero. Based on this as- 
sumption, we plot (as in the example shown in Figure 5) the 
cross-correlation of bs (vi,V2) (along the y-axis) against the as- 
sociated product of mean intensities 5 ave (vi) 5 ave (v2). From 
this plot an estimate of of, or the slope of the expected lin- 
ear dependence, can be derived. The velocity channels where 
the signal-to-noise ratio of the mean intensity S ave (v) is less 
than 3 are excluded from this analysis. Naturally, the num- 
ber of useful channels (M, typically about 5-6) based on this 
criterion is significantly reduced, but the resultant number of 
correlation pairs [M(M- l)/2] is large enough. The observed 
scatter about the linear dependence arises from the terms with 
f v and/or n, resulting in the uncertainty in the estimation of its 
slope. The two straight lines about the best-fit line, all passing 
through the origin, indicate the ±<r uncertainty in the slope. 
With of estimated in this manner, the "broadband" modula- 
tion contribution can now be removed from the observed vari- 
ance. 

In practice, any intrinsic fluctuation which may other- 
wise be uncorrelated between adjacent channels may con- 
tribute to some (positive) correlation due to the finite veloc- 



ity resolution of the measured spectra (0.11 km s -1 at 1665 
MHz). Since such contributions will mimic correlations due 
to "broadband" modulation, these contributions can lead to an 
over-estimation of of, and as a result, the contribution from 
"narrowband" fluctuations will be correspondingly underesti- 
mated. Any intrinsic fluctuation with a correlation bandwidth 
wider than the channel width will also be underestimated. If 
these effects are in fact significant, our estimates of the modu- 
lations index (or variance) associated with intrinsic variability 
can then be viewed as lower limits. As one possible measure 
against over-estimation of of, we limit its value, if needed, 
such that of (v)) is never negative (see Equation 5). 

To assess the possible contamination due to these effects, as 
well as the robustness of the correlation procedure, we have 
repeated the analysis after normalization of each profile in 
the observed dynamic spectra with respect to (1) the velocity- 
averaged intensity, and (2) the intensity in a reference channel 
corresponding to the peak of the average profile. The esti- 
mates of <j c in these two cases are close to zero, suggesting 
that the above discussed contaminations are not significant. 
Also, the estimates of the intrinsic modulation are found to be 
consistent across the three methods, indicating that the esti- 
mate is not sensitive to whether or how the "broadband" vari- 
ation was modified before applying the correlation procedure. 

We have also examined the performance of the three meth- 
ods when used separately. Although, they provide a consis- 
tent accounting (and removal) of the "broadband" modulation 
contribution, there are significant differences. Normalization 
with respect to the peak channel intensity is based on an im- 
plicit a priori expectation that the intrinsic variability in that 
reference channel is absent, and any observed variability is 
therefore only of "broadband" nature. Invalidity of this ex- 
pectation can lead to a systematic over-estimation of the vari- 
ance in the other channels, depending on the modulation index 
of the "narrowband" variation (including system noise) in the 
reference channel. However, for the channels adjacent to the 
reference channel sharing any common "narrowband" vari- 
ability, the resultant variance would be underestimated. On 
the other hand, normalization based on the velocity-average 
intensity makes no assumption about absence of intrinsic vari- 
ability in any particular channel, but does implicitly assume 
that the line-integral is intrinsically constant. Such an as- 
sumption has no physical justification. Hence, any "narrow- 
band" variability, apart from its magnitude being underesti- 
mated, contaminates other channels with anti-correlated vari- 
ations and a consequent increase in the variance. In both of 
these approaches, the corrected dynamic spectra are available, 
and can be examined using fluctuation spectral analysis for 
estimating the temporal scales of the remaining variability. In 
contrast, the correlation-based estimation and removal of the 
"broadband" contribution to the variance only produces time- 
averaged quantities, and hence any further temporal/spectral 
analysis can not be attempted. However, the correlation-based 
approach is the most unbiased in comparison with the two 
methods based on normalization. Therefore, for a determina- 
tion of the variance associated with the "broad- and narrow- 
band" variability, we have used the correlation-based proce- 
dure described above, and have only employed the normal- 
ization by intensity in the reference peak channel for the pur- 
poses of temporal/spectral analysis (e.g. ^6). 

The results from our variance analysis on all of the 6 data 
sets (consisting of dynamic spectra for RR & LL polariza- 
tion for the sources in Table 1) are summarized in Figure 
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Product of mean intensities (Jy/beom) 2 

FIG. 6.— An example, based on W30H 1665-MHz line (RR) data, to il- 
lustrate our procedure to estimate the possible contribution from "broadband" 
modulation. Here, we examine the cross-correlation between variations ob- 
served in every pair of channels versus the product of mean intensities cor- 
responding to the respective channel pair. The observed trend is shown by 
the dashed line. Only those channels with their mean intensities having S/N 
greater than 3 are included. The slope of the expected linear dependence pro- 
vides an estimate of (see Equation 6), where a c is the modulation index 
associated with the "broadband" modulation that is shared by all channels. 
The observed scatter about the linear dependence arises in general from any 
"narrowband" variations including the measurement noise. The middle of the 
three solid lines, always passing through the origin (0,0), indicates the best fit 
to the data, while the two other lines correspond to the ±1 — a uncertainty in 
the slope. 



showing in each case the profiles of average line inten- 
sity (S ave ), as well as the standard deviations associated with 
the observed variability (<7 & s (v)), possible "broadband" mod- 
ulation (<Xb(v) = &cSave(v)) an d the measurement uncertainties 
(cr n ). As already described and illustrated in Figure [6] the 
"broadband" modulation index a c is estimated based on the 
cross-correlation analysis. The quantity cr^ b — a b , where the 
"broadband" modulation contribution is removed, represents 
the observed "narrowband" variance, which includes the nom- 
inally expected contribution o\ from the measurement noise. 
Also, we examine the ratio of the observed to the expected 
variances, R = (o^ bs — cP^)ja\, for any significant deviations 
from the expected value of unity. In other words, any statis- 
tically significant residual variance, i.e. af = a^ bs — of— a\, 
must then be "narrowband" in nature, and thus may be as- 
sumed to be intrinsic to the source. As another measure, we 
also compute the modulation index associated with the pos- 
sible intrinsic "narrowband" variability, a v = Ui/S ave , where 
the signal-to-noise ratio of S ave estimate is three or more. The 
profiles of the ratio R and a v are also displayed in Figure0 

A number of salient facts are evident from FigureQ Firstly, 
statistically significant variations are indeed detected, even af- 
ter accounting for the "broadband" modulation and measure- 
ment uncertainties. These variations, reflected by the excess 
variance at certain velocities, are "narrowband" in nature, and 
hence intrinsic to the source. Moreover, these variations are 
apparent at velocities away from the line peak, rather than at 
the peak. Within the statistical errors, the peaks do not appear 
to exhibit significant variations. The variations observed in 



the shoulders are detectable with signal-to-noise ratio of three 
or more. 

6. VELOCITY-RESOLVED FLUCTUATION SPECTRA, AND 
TIME-SCALES OF INTRINSIC VARIABILITY 

The basic data used here also are the dynamic spectra span- 
ning 12.5 hour and obtained as chronologically arranged line 
profiles from one minute integrations. The aim of the fluctua- 
tion spectral analysis described below is twofold. The first is 
to examine the nature of the fluctuation spectra. Then, given 
the fluctuation spectra, to estimate the timescale of flux den- 
sity variations arising from "narrowband" intrinsic variations. 
Thus it is essential to eliminate contribution from any "broad- 
band" variations. To ensure this, each of the spectral profiles 
in the dynamic spectra is normalized with respective intensi- 
ties in the reference channel defined by the peak in the aver- 
age line profile. The justification for this normalization comes 
from the variance analysis where we concluded that the 
peak of the line profile does not exhibit any statistically sig- 
nificant "narrowband" variation. We can also rule out any 
overestimation of variability in other channels due to the nor- 
malization. 

These "corrected" dynamic spectra for each of the sources 
form the input data for the velocity-resolved fluctuation spec- 
tral analysis. This analysis is very similar to the longitude- 
resolved fluctuation spectral analysis, a well known tool in 
pulsar emission studies (Backer 1973, Deshpande & Rankin 
1999). Here, we compute fluctuation power spectrum for each 
of the velocity channels separately, where the single tem- 
poral sequence of intensities in a given channel is Fourier 
transformed and the power at each of the fluctuation frequen- 
cies computed. The results corresponding to all the velocity 
channels for a given line source are displayed together as a 
two-dimensional display of fluctuation power as a function 
of velocity and fluctuation frequency, along with a velocity- 
averaged power spectrum (see Figure[8}. 

Any statistically significant power observed in these fluc- 
tuation spectra can be interpreted as due to variability that is 
necessarily "narrowband", and hence, intrinsic to the source. 
The following aspects are clearly evident from the spectra in 
Figure[8] The fluctuation power is generally higher in the line 
channels compared to those well away from the line emis- 
sion, as would be expected from the measurement noise that 
will be proportional to the system temperature including the 
line intensity. Moreover, another thing that is apparent is 
that the fluctuation power toward lower frequencies is more 
than that at higher fluctuation frequencies. We examine also 
the spectrum computed by averaging the fluctuation spectra 
across the velocity channels. With the equivalent line width 
of only a couple of channels, given that the contribution of the 
line channels overwhelms in this average, the resultant spec- 
trum has benefited only correspondingly from the averaging 
across velocity. Hence the ratio of the average power to the 
uncertainties (i.e. the signal-to-noise ratio) at any fluctuation 
frequency is rather small (~1.4), except when the fluctuation 
power in the line channels is small (e.g. as seen at the high 
fluctuation frequencies in some of the average spectra). Hence 
we do not consider the apparent fine spectral structure as sig- 
nificant, but rather examine and estimate the smooth trends 
across the spectrum, since they will have a much improved 
significance in accordance with the smoothing scale. The 
fluctuation power level at higher frequencies is consistently 
low, and corresponds largely to the measurement-noise contri- 
bution that is generally expected to be "white" (or uniform) in 
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FIG. 7. — Results of the variance analysis of the dynamic spectral data on the sources listed in Table 1. The plots corresponding to the LL & RR polarization 
channels are shown on the left and the right columns, respectively, The top, middle and the bottom pairs show the results for the 1612, 1665 and 1720 MHz 
transitions, respectively. In each of the six plots, the top panel shows the average intensity profile (solid line), along with the profiles of standard deviations 
associated with the observed variance (dashed line), measurement noise (dotted line) and the "broadband" modulation (dot-dashed line), respectively. The latter 
is a scaled version of the average intensity profile, where the scale factor a c is estimated using the cross-correlation procedure in ^5] For clarity, all three profiles 
for the standard deviations are amplified by a factor of ten in this display. The middle panels display profiles of the ratio (R) of the observed "narrowband" 
variance to its expected value. The observed "narrowband" variance is simply the observed variance minus that associated with the "broadband" modulation, 
while the expected "narrowband" variance is that from the measurement noise The three curves correspond to the mean value of the ratio R, and with 

±1 — <r deviations from this mean. The nominally expected values of R is unity within the indicated uncertainties, and any significant excess is interpreted as 
due to intrinsic "narrowband" variability. The modulation index (MI) associated with the latter, along with its ±1 — a bounds are given in the bottom panel. The 
modulation index (MI) is computed as the ratio of the RMS fluctuations attributable to the intrinsic variability to the mean intensity Save, when both of them are 
significantly greater than zero. All of the estimates presented above are based on the data spanning 12.5 hours. 
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its spectral character. Thus the overall increase in the fluctua- 
tion power toward the lower frequency portion is significant in 
most cases (except for the 1665 line), where the power is typi- 
cally about two to three times that at the higher frequency end 
of the average spectrum. The latter defines the reference noise 
floor, across which the RMS variation may be estimated and 
used to assess the significance of the contrast of fluctuation 
power levels between the two frequency ranges. For example, 
the increase in the average power toward lower frequencies 
for the 1612-MHz RR line is about seven times the RMS de- 
viation in the power at the higher frequencies. This factor is 
somewhat lower for other lines/polarization channels, and is 
close to zero for the 1665 RR line. 

We observe significant relative fluctuations up to a fluc- 
tuation frequency of ~ 10~ 3 Hz, which corresponds to time 
scales of > 15-20 min. Given that the normalization proce- 
dure has removed any extrinsic variability that is expected to 
be "broadband", we associate this observed time scale with 
the intrinsic variability of the OH maser lines. 

7. DISCUSSION AND SUMMARY 

In this work, we have conclusively demonstrated the pres- 
ence of "narrowband" intrinsic variations in these W3(OH) 
maser sources. These variations seem to have a typical time 
scale of about 15-20 minutes or longer, indicating that the 
faster time scale may correspond to the light travel time of 
the maser pumping column. 

A very important aspect, which is worth stating again, is 
that the observed intrinsic variations that are clearly sepa- 
rable from other contaminants are necessarily "narrowband" 
in nature. All the other variations such as the instrumental 
gain variations and the interstellar diffractive scintillation, al- 
though their time scales may be of the order of a few tens 
of minutes, are distinguishable from the intrinsic variabilities, 
mainly because of their broad correlation bandwidth. The ex- 
pected differential fluctuation within the line profile due to 
interstellar scintillation is only of the order of ^0.5%, much 
smaller than the observed intrinsic variability, whose magni- 
tude was as much as 100% (modulation index), in the "tail" 
portions of the spectral lines. 

There is one source of systematic narrowband error that 
could have potentially influenced our conclusions. In our 
VLBI measurements, if the Earth's motion were not correctly 
compensated for, certain systematic but spurious temporal 
variations in the line profiles would be induced due to dif- 
ferential Doppler shift. However, we would then expect the 
variations seen on the two sides of the line profile peak to 
be anti-correlated. We have examined carefully the relevant 
cross-correlations, and do not find any statistically significant 
signature such an anti-correlation, clearly indicating that our 
data are free of such artefact. 

An important aspect of these maser sources is that their typ- 
ical measured angular diameter is ^3 mas. With the distance 
to the star forming region of 2 kpc, this corresponds to a trans- 
verse distance scale of ~ 6 AU. As we have seen from Figure 
|8j the fluctuation frequencies seen in the spectra are < 10~ 3 Hz 
(fluctuation time scale, 7) ~ 1000 s or longer). If indeed this 
time scale reflects the dimension of the source based on light 
travel time arguments, then the implied longitudinal spatial 
scale of the maser column would be 2-3 AU. It is important 
to note that the apparent transverse spatial dimension of the 
source as measured from the angular size is comparable to its 
longitudinal dimension implied by the fluctuation timescale, 
even though there seems to be no a priori basis for the com- 



parison, let alone the agreement. However, if indeed these 
orthogonal dimensions of the source are expected to compa- 
rable, it would imply that the VLBA observations may have 
actually resolved the intrinsic source size of the OH masers in 
W3(OH). Then the possible contribution of angular broaden- 
ing caused by interstellar scattering to the apparent size of the 
source is minimal. This also implies that the expected decor- 
relation bandwidth of interstellar scintillation may be much 
wider than 100 kHz, further reinforcing the validity of our as- 
sumptions. However, if the scatterer is closer to the source 
instead of the midway location which is implicit in equation 
1, the decorrelation bandwidth would be correspondingly nar- 
rower. We assess whether that may be the case, based on the 
two important estimates we have at hand, namely the mea- 
sured upper-limit for scatter broadening Oh and the observed 
time scales (say, 7^) of the narrow-band variability, if due 
to ISM, since that is our concern. How far are we from the 
naive assumption of uniform scatterer (or a strong scatterer 
mid-way) will be reflected in both of these parameters, which 
depend differently on the relative location of the scatterer (see 
for example, eq. 3 and 8 of Deshpande & Ramachandran 
1998). If D s & D are the distances of the scatterer and of the 
observer, respectively, as measured from the source, it is easy 
estimate an upper limit to (D/D s -1) using relevant expres- 
sions, and assuming that the observed time-scales are 1000 s 
or longer, & a wavelength of 18 cm. We find that (D/D s -1) 
should be less than or equal to 20/(6*// V, sm ), where the scatter 
broadening Oh is in mas, and the ISM velocity V,„„ is in km 
s" 1 . A value of 3 mas for Oh implies (D/D s - 1) < 6.6/V/ jm . 
Noting that V, sm will have the same lever-arm factor as source 
velocity does, and considering typical values of V,„„, we con- 
clude that (D/D s - 1 ) is close to unity, if not smaller. Hence, 
the decorrelation bandwidth that we estimate naively is most 
likely an under-estimate, given the above and that Oh may al- 
ready be an over-estimate of the scatter broadening. 

As already mentioned in the introduction, variability of as- 
trophysical masers on longer time scales — weeks to months 
— observed before by several investigators, may be due to 
changes in the number density of relevant molecules or phys- 
ical conditions in the maser column. On the other hand, very 
short time scale rapid variations have also been seen. Salem 
& Middleton (1978) suggest a model in which a sudden on- 
set of a pumping mechanism can cause rapid quasi-periodic 
fluctuations in the observed flux density, and predict fluctua- 
tions with time scales of a day or so, with ~25% modulation 
index. These fluctuations may either correspond to propaga- 
tion of a radiative or a collisional disturbance. In the former 
case the disturbance travels with the speed of light, and in 
the latter with a typical speed of ^10 km s" 1 or so. For ra- 
diative propagation of disturbance, with the typical flux den- 
sity of the observed lines, the brightness temperature comes 
to ~ 10 10 K. For a collisional disturbance, the source dimen- 
sion is only 7} ■ x 10 km s" 1 « 10 7 meters (angular diameter 
of ^0.1 /larcsec). This corresponds to a brightness tempera- 
ture of ~ 10 20 K (see also Clegg & Cordes 1991). In the lat- 
ter case, since the implied intrinsic angular diameter is only 
~ 1 /iarcsec, the observed angular diameter is predominantly 
due to interstellar scattering. For an object at 2.2 kpc (dis- 
tance of W3(OH) complex), a scatter broadening of 2-3 arc- 
sec seems to be very large at the 18-cm wavelength. VLBI 
observations of several nearby (distance less than 2-3 kpc) 
pulsars at 327 MHz have shown a broadening of 10-20 mas, 
with an exception of the Vela pulsar, PSR B0833^-5 (Gwinn, 
Bartel & Cordes 1993). In the case of Vela, the excess an- 
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FIG. 8. — Results of the Velocity-resolved fluctuation spectral analysis. The pair of plots from top to bottom correspond to the 1612, 1665 and 1720 MHz 
lines, respectively. In each of these pairs, the left and the right plots are for the LL, and RR polarization channels, respectively. Within each plot, the left side 
panels shows the average line profile (flux density in units of Jansky/beam) as a function of velocity obtained from the full observation. The central panel shows 
the velocity resolved fluctuation spectra, where each horizontal row is the power spectrum of the time series from the corresponding velocity channel. Each time 
series is normalized with respect to the time series for the reference peak channel. Naturally, this results in a null spectrum for the reference channels. As the color 
bar shows, red corresponds to the maximum power, and white corresponds to zero power. The bottom panel gives the velocity-averaged fluctuation spectrum. 
In most cases, significant fluctuation power is apparent at lower frequencies in comparison with the noise floor at the higher frequency end of the spectra. As 
discussed in the main text, the observed "narrowband" variability on time scales of 15-20 minutes or longer is to be understood as intrinsic to the maser sources. 
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gular broadening is due to the presence of strong scattering 
screen (Deshpande & Ramachandran 1998). Assuming a Kol- 
mogorov density irregularity spectrum, with the wavelength 
dependence of A 2,2 , the expected angular broadening is only 
< 0.5 mas at 1665 MHz. Based on these considerations, we 
argue against any significant over-estimation of the source 
size of W3(OH) due to scatter broadening. 

It has been suggested by Eli tzur (1991) that short time scale 
intensity variations can be produced by variations in the maser 
level population at small length scales. These fluctuations 
produced in the unsaturated medium of the maser core give 
rise to the spectrum of flux density variability observed. In 
our case, with the length scale of 2-3 AU (derived based on 
the light travel time arguments) represents an average "seed" 
length for such a fluctuation. This suggestion may well be 
the reason for the short time scale fluctuations that we have 
observed. 

As mentioned earlier, the longitudinal dimension 2-3 AU 
estimated from the observed time scales of intrinsic variabil- 
ity compares well with the transverse spatial dimension of 6 
AU estimated assuming that the apparent angular size to be 
approximately the intrinsic size of the source. However, for 
estimating brightness temperature the source we will use only 
its transverse size. A typical observed flux density of some 
10 Jy at the line peaks (e.g. for 1612/1720 lines), implies a 
peak brightness temperature of ~ 2 x 10 13 K. A five times 
higher value, i.e. ~ 10 14 K, is implied by the correspondingly 
brighter peaks of the 1665 lines. 

The intrinsic intensity variations that we observe are partic- 
ularly confined to the "shoulders" of the lines and away from 
the line peak (see Figure[7}. The peak of the line profiles does 
not exhibit any statistically significant variation, in all of the 
six data sets. This behavior is not at all surprising if the peak 
of these line emissions corresponds to "saturated" maser ac- 
tion, and the rest of the line profile exhibits unsaturated emis- 
sion. In any case, the absence of any intrinsic variability at the 
line peak enables use of the corresponding channel intensities 
to calibrate out any "broadband" variations from the dynamic 
spectrum. 

Another interesting aspect apparent from the results in Fig- 
ureQis that the statistically significant narrowband variations 
that we observe at the raising and falling edges of the line pro- 



files are not identical between the Zeeman pairs (RR and LL 
components). For instance, the difference between the mod- 
ulation indices of the two Zeeman components of 1720 MHz 
lines is clearly seen to be several times the root mean square 
noise level. The physical reason for this behavior needs to be 
explored. Given the birefringent nature of the medium close 
to or at the origin of the radiation, as well as interstellar propa- 
gation medium, effects of refraction on the apparent visibility 
of the two Zeeman components remain to be explored. 
To summarize, we list our main conclusions as follows. 

• The combination of our cross-correlation procedure and 
the variance analysis provides an effective tool for esti- 
mation and elimination of the "broadband" contribution 
from instrumental effects and the interstellar diffractive 
scintillation, and thus for identifying variations that are 
intrinsic and "narrowband" in nature. 

• Significant intrinsic "narrowband" variability is ob- 
served over most regions of the line profiles except at 
and about the line peak, suggesting "saturated" maser 
emission at the peak. 

• The velocity-resolved fluctuation spectra reveal that the 
time scales of the significant intrinsic variability are 15- 
20 minutes or longer. 

• Based on light travel time argument, the intrinsic vari- 
ability time scale implies a longitudinal spatial scale of 
the maser column to be about 2-3 AU. 

• The apparent angular sizes of the sources are unlikely to 
be significantly affected by interstellar scattering. The 
implied transverse size of the source is comparable to 
its longitudinal dimension. 

• The peak brightness temperatures for the maser sources 
range between ~2x 10 13 and ~ 10 14 K. 
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